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The PHENIX 0 and STAR |Q, ^ Collaborations at the Relativistic Heavy Ion Collider 
(RHIC) have recently presented data for production of heavy quarks in pp and d+Au collisions 
at s/Snn = 200 GeV, both in the form of explicitly reconstructed D mesons, and as an electron 
spectrum from semi-leptonic decay of the heavy hadrons. These data must be compared to QCD 
predictions, to establish the extent to which they can be successfully described, before moving 
forward and trying to determine the presence of dense matter effects in nucleus-nucleus collisions. 

Recent improvements in heavy quark production theory and experimental measurements at 
colliders, especially for bottom production, have shown that the perturbative QCD framework 
works rather well, see Refs. [Q], provided a number of precautions are taken in performing the 
phenomenological analysis and the proper modem tools are employed. In particular, the advan¬ 
tages of comparing the theoretical prediction and the data directly at the level of the experimental 
observables have been clearly outlined. 

The purpose of this work is to repeat these analysis for RHIC, and to provide therefore a 
benchmark prediction against which the data can be compared. 

We calculate the transverse momentum (px) distributions of charm and bottom quarks, the 
charm and bottom hadron distributions resulting from fragmentation and, finally, the electrons 
produced in semi-leptonic decays of the hadrons |^. Theoretical uncertainties, estimated as exten¬ 
sively as possible, constitute an intrinsic component of the prediction. Our final result is thus not a 
single curve, but rather an uncertainty band which has a reasonably large probability of containing 
the ‘true’ theoretical prediction. 

The theoretical prediction of the electron spectrum includes three main components: the px 
and rapidity distributions of the heavy quark Q in pp collisions at s/S = 200 GeV, calculated 
in perturbative QCD; fragmentation of the heavy quarks into heavy hadrons, Hq, described by 
phenomenological input extracted from data; and the decay of Hq into electrons according to 
spectra available from other measurements. This cross section is schematically written as 
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where the symbol (g) denotes a generic convolution. 

The distribution Ed^o{Q)/dp^Q is evaluated at Fixed-Order plus Next-to-Leading-Log (FONLL) 
level, implemented in Ref. 0. In addition to including the full fixed-order NFO resulf 0, fhe 
FONFF calculafion also resums 0 large perfurbafive ferms proportional fo a”log^(px/m) fo all 
orders wifh nexf-fo-leading logarifhmic (NFF) accuracy (i.e. k = n,n — 1) where m is fhe heavy 
quark mass. The perfurbafive paramefers are m and fhe value of fhe sfrong coupling, a^- We fake 
mc = 1.5 GeV and mi, = 4.75 GeV as cenfral values and vary fhe masses in fhe range 1.3 < me < 1.7 
GeV for charm and 4.5 <mb <5 GeV for boffom fo esfimafe fhe mass uncerfainfies. The five-flavor 
QCD scale is fhe CTEQ6M value, A^^^ = 0.226 GeV. The perfurbafive calculafion also depends on 
fhe faclorizafion (pf) and renormalization (Pr) scales. The scale sensifivify is a measure of fhe 
perfurbafive uncerfainfy. We fake p^ f = po = s/p^ + m^ as fhe cenfral value and vary fhe fwo 
scales independenfly wifhin a ‘fiducial’ region defined by p^ p = ^RpPo wifh 0.5 < ^Rp < 2 and 
0.5 < ^r/^f < 2 so fhaf = {(1,1), (2,2), (0.5,0.5), (1,0.5), (2,1), (0.5,1), (1,2)}. The 

envelope confaining fhe resulfing curves defines fhe uncerfainfy. The mass and scale uncerfainfies 
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Figure 1: Left-hand side: The theoretical uncertainty bands for c quark and D meson pj distributions in 
pp collisions at s/S — 200 GeV, using BR(c ^ D) = 1. The final and preliminary [|^] STAR dn-Au data 
(scaled to pp using Abin = 7.5) are also shown. Right-hand side: The same for b quarks and B mesons. 


are then added in quadrature. These inputs lead to a FONLL total cc cross section in pp collisions 
of = 256^j 4 g ph at s/S = 200 GeV, and to a total bb cross section of 1.87+Qgy ph. 

The fragmentation functions, D{c D) and D{b —> B), where D and B indicate a generic 
admixture of charm and bottom hadrons, are consistently extracted from data in the context of 
FONLL [^[T^. Using the Peterson et al. fragmentation function, with standard parameter choices 
Ec — 0.06 ± 0.03 and Sf, ^ 0.006 ± 0.003, does not provide a valid description of fragmentation in 
FONLL, since the hadronization is too soft. In fact. Ref. showed that replacing the Peterson 
fragmentation description with an appropriate one constitutes one of the main improvements which 
help reconcile the bottom transverse momentum distribution measured at the Tevatron with the 
theoretical prediction. 

The measured spectra for primary B e and D e decays are modeled and assumed to be 
equal for all bottom and charm hadrons, respectively. The contribution of electrons from secondary 
B decays, B ^ D ^ e, was also included but is mostly negligible. The decay spectra are normalized 
using the branching ratios for bottom and charm hadron mixtures: BR(B e) = 10.86 ±0.35%, 
BR(D -^e) = 10.3 ± 1.2%, and BR(B ^ D ^ e) = 9.6 ± 0.6%. 

The left-hand side of Fig. |T] shows the theoretical uncertainty bands for c quarks and D mesons, 
obtained by summing the mass and scale uncertainties in quadrature. The band is broader at low 
Px due to the large value of and the behavior of the CTEQ6M parton densities at low scales as 
well as the increased sensitivity of the cross section to the charm quark mass. 

Figure § compares the RHIC data to the total uncertainty band for D ^ e, B —> e and B 
D —> e decays to electrons. The upper and lower limits of the band are obtained by summing the 
upper and lower limits for each component. It is worth noting that, while for the central parameter 
sets, the B —> e decays begin to dominate the D —> e decays at ~ 4 GeV, a comparison of the 
individual bands (not shown) shows that the crossover may occur over a rather broad range of 
electron pj. The relative c and b decay contributions may play an important part in understanding 
the electron Raa in nucleus-nucleus collisions which seems to suggest strong energy loss effects on 
heavy flavors Q,0. 

In conclusion, we have performed a phenomenological analysis of heavy quark production in 
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Figure 2: The final prediction for the theoretical uncertainty band of the electron spectrum from charm and 
bottom in pp collisions. Data from PHENIX [|^ and STAR (final and preliminary [^) are also shown. 


s/S = 200 GeV pp collisions at RHIC. The results are presented in the form of a theoretical uncer¬ 
tainty band for the transverse momentum distribution of either bare charm (bottom), D (B) mesons, 
or electrons originating from the decay of charm and bottom hadrons. This band stems, at the per¬ 
turbative level, from a next-to-leading order result improved by next-to-leading log resummation at 
large transverse momenta. These result should not be multiplied by any K factor before comparison 
with data. Rather, agreement within the uncertainties of the measurements will support the applica¬ 
bility of standard QCD calculations to heavy quark production at RHIC. Alternatively, a significant 
disagreement will suggest the need to complement this evaluation with further ingredients. 
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